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■ Abstract 

Existing predictions for the branching ratio for — > KK'y via <p — > £7 (where 
^ ■ S denotes one of the scalar mesons / (975) and a (980)) vary by several orders of 

magnitude. Given the importance of these processes for both hadron spectroscopy 
and CP-violation studies at cf) factories (where <ft K°K 0/ -f poses a possible 
background problem), this state of affairs is very undesirable. We show that 
the variety of predictions is due in part to errors and in part to differences in 
modelling. The latter variation leads us to argue that the radiative decays of 
these scalar states are interesting in their own right and may offer unique insights 
into the nature of the scalar mesons. As a byproduct we find that the branching 
ratio for <p — > K°K 0/ -f is < 0(10~ 7 ) and will pose no significant background to 
proposed studies of CP-violation. 
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1 Introduction 



There are predictions in the existing literature for the branching ratio for 
— > KK'y via — > £7 (where 5 denotes one of the scalar mesons S* (now called 
/o(975)) or 5 (now called a (980)) that vary by several orders of magnitude [1- 
5]. Clearly not all of these predictions can be correct! Given the importance of 
these processes for both hadron spectroscopy and CP-violation studies, this state 
of affairs is clearly undesirable. Moreover, in view of the impending factory, 
DA$NE [6], and other developing programmes [7], there is an urgent need to 
clarify the theoretical situation. 

The scalar mesons (i.e., mesons with J PCn = ++ ) have been a persistent 
problem in hadron spectroscopy^ We shall show in this paper that the radiative 
decays of the meson to these states can discriminate among various models of 
their structure. In addition to the spectroscopic issues surrounding the scalar 
mesons, there is a significant concern that the decay K°K 0, -f poses a possible 
background problem to tests of CP-violation at future factories: the radiated 
photon forces the K°K° system to be in a C-even state, as opposed to the C- 
odd decay — > K°K°. Looking for CP-violating decays in — ► K°K° has been 
proposed as a good way to measure e' /e [10], but because this means looking for 
a small effect, any appreciable rate for — ► K°K 0/ -f (namely, a branching ratio 
— > K°K 0, y > 10~ 6 ) will limit the precision of such an experiment. Estimates 
[4] of the non-resonant — > K°K 0/ -f rate give, in the absence of any resonant 
contribution, a branching ratio of the order of 10~ 9 , far too small to pose a prob- 
lem. The uncertainty in the theoretical estimates, and the potential experimental 
ramifications, arise due to the presence of the scalar mesons /o(975) and a (980), 
which are strongly coupled to the KK system. Estimated rates for the resonant 
decay chain — > S + j, followed by the decay S — > K°K°, vary by three orders of 
magnitude, from a branching ratio of the order of 10 -6 down to 10~ 9 . These vari- 
ations in fact reflect the uncertainties in the literature for the expected branching 
ratio for — ► 5*7 which vary from 10~ 3 to 10 -6 [11] . Here we concentrate on this 
resonant process. 

We shall show that the variability of the predictions for — > 5*7 is due in 
part to errors and in part to differences in modelling. On the basis of this model 
dependence, we argue that the study of these scalar states in — > S7 may offer 
unique insights into the nature of the scalar mesons. These insights should help 

*For an historical perspective see Ref. [8]; for a more recent study see Ref. [9]. 
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lead in the future to a better understanding of not only quarkonium but also 
glueball spectroscopy. As a byproduct we predict that the branching ratio for 
<f) -> K°K°-f is < 0(l(r 7 ) (i.e., the branching ratio for <p -> S 7 is < 0(10~ 4 )) 
and will pose no significant background to studies of CP-violation at DA$NE. 

2 Probing the nature of the scalar mesons below 1 GeV 



The scalar mesons are spectroscopically interesting for several reasons. One 
is that, while agreeing on little else, it is an essentially universal prediction of 
theory (lattices, bags, flux tube models, QCD sum rules, . . . ) that the lowest- 
lying glueball has scalar quantum numbers and a mass in the 1.0 - 1.5 GeV mass 
range. Clarifying the presently confused nature of the known ++ mesons may be 
pivotal in the quest to identify this glueball. Another is the possibility that the 
two best known [12] scalar mesons, the /o(975) and the ao(980), are qqqq states. 
The original proposal [13] for this interpretation, based on the bag model, also 
predicted many other states which have not been seen (although this shortcoming 
is now understood to some degree [14]). The qqqq interpretation of these two 
states was later revived in a different guise within the quark potential model as 
the U KK molecule" interpretation [15]. Since providing a test of this particular 
interpretation is one of the main results to be presented here, we first briefly 
elaborate on these two models of multiquark states. 

In the naive bag model the qqqq states consist of four quarks confined in a 
single spherical bag interacting via one gluon exchange. It is obvious that such a 
construction will lead to a rich spectroscopy of states. Although it is not clear how 
to treat or interpret the problem of the stability of this spectrum under fission into 
two bags [14], it is very interesting that the dynamics of this model predicts that 
the lowest-lying such states will (in the SU(3) limit) form an apparently ordinary 
( "cryptoexotic" ) nonet of scalar mesons. It is, moreover, probable that a better 
understanding of bag stability could solve both the problem of unwanted extra 
predicted states and also a problem with the do itself: in the naive model it can 
"fall apart" into tti] so that it is difficult to understand its narrow width, given 
the presently accepted pseudoscalar meson mixing angle (see footnote 22 in the 
first of Refs. [13]). In the absence of an understanding of how to overcome these 
difficulties, we will not consider the bag picture further in this paperQ. 

^See, however, Refs. [16] for a possible way out of the a — > ttij problem. 
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In the potential model treatment [15] it is found that the low-lying qqqq sector 
is most conveniently viewed as consisting of weakly interacting ordinary mesons: 
the resulting spectrum is normally a (distorted) two particle continuum. Within 
the ground state u,d,s meson-meson systems, the one plausible exception to this 
rule is found in the KK sector (i.e., the KK channel and those other channels 
strongly coupled to it): the L = (i.e., J PCn = ++ ) spectrum seems to have 
sufficient attraction to produce weakly bound states in both 1 = and 1=1. 
There are a number of phenomenological advantages to the identification of these 
two states with the / (975) and a (980). Among them are: 

1) It is immediately obvious why the / (975) and a (980) are found just below 
KK threshold: they bear much the same relationship to it that the deuteron bears 
to np threshold. 

2) The problem of the / and a widths is solved. If these states were 3 Pq 
quarkonia with flavours corresponding to uo and p (as suggested by their degener- 
acy), then r(/ — > 7T7r)/r(a — > nr]) would be about 4 in contrast to the observed 
value of about \. At least as serious is the problem in the quarkonium picture 
with the absolute widths of these states: models [17-19] predict, for example, 

r(/o^7T7r) ~ (3-6)^^(^)5) (2.1) 
~ 500 - lOOOMeV (2.2) 

versus the observed partial width of 25 MeV. We have already noted the problem in 
the bag model qqqq interpretation with a Q — > nrj. In the KK molecule picture the 
narrow observed widths are a natural consequence of weak binding: (KK)i =0 — > 
7T7T and (KK) I= i — > nr] occur slowly because the KK wavefunction is diffuse. 

3) Both the / and a seem to bear a special relationship to ss pairs: their 
KK "couplings" are very large and they are observed in channels which would 
violate the Okubo-Zweig-Iizuka (OZI) rule [20] for an cu, p -like pair of states [21]. 

4) The 77 couplings of the fo and ao are about an order of magnitude smaller 
than expected for 3 P quarkonia [22], but consistent with the expectations for KK 
molecules [23]. 

Although these observations argue against the viability of the 3 Po quarkonium 
interpretation of the / (975) (and probably also the a (980)), they are insufficient 
to rule it out completely. (Moreover, a unitarized variant of the quark model [24], 
in which the scalar mesons are strongly mixed with the meson-meson continuum, 
avoids several of these problems. In addition to this conservative alternative, the 
recent analysis of Ref. [9] has raised the possibility that the / (975) is really a 
combination of two effects, one of which is a candidate for a scalar glueball.) 
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The main purpose of this paper is to point out a simple (and to us unexpected) 
experimental test which sharply distinguishes among these alternative explana- 
tions. We show that the rates for — > fo (975)7 ~~ n7T l an d — > ao (980)7 ~~ 7r? 77 
in the quarkonium, glueball, and KK molecule interpretations differ significantly; 
furthermore, the ratio of branching ratios 

-> a (980)7 
-+ /o(975) 7 

also may prove to be an important datum in that it can have a model-dependent 
value anywhere from zero to infinity (see Table 2)! 

In the quarkonium interpretation, — > f (975)7 an d — * a (980)7 are s i m_ 
pie electric dipole transitions quite similar in character to several other mea- 
sured electric multipole transitions, including not only the light quark transitions 
a 2 (1320) -> 7T7, AT* (1420) -> A7, ai(1275) -> vr 7 , and 6i(1235) -> 7r 7 , but also 
such decays as Xco —> "07 an d Xbo — * Y7. From the comparison between theory 
and experiment given in Ref. [17], we expect that the quark model predictions 
for these processes given in Table 1 are reliable to within a factor of two. Thus if 
the /o is an ss quarkonium, the branching ratio for — > S7 would typically be of 
the order of 10~ 5 . 

If the /o(975) is a glueball (in Ref. [9] there is a glueball component of the 
U S* effect", dubbed the 5*1(991), which couples to mv and is responsible for the 
resonant behaviour seen in hit phase shift analyses; the other component, dubbed 
the 52(998), is practically uncoupled to nn) then one would naturally expect 
— > fo (975)7 ~" *■ 7r7r 7 to be even smaller than in the quarkonium interpretation 
since the decay would be OZI- violating. The remarks made above on the strong 
decay widths of the quarkonium states would suggest that quarkonium - glueball 
mixing, through which we presume the OZI-violation would proceed, must be 
small for the / (975) to remain narrow. Thus we can crudely estimate the glueball 
- quarkonium mixing angle to be less than [r(/ — ► kit)/T( 3 P — > itit)]^ so that if 
fo (975) is a glueball 

r(0 - f (glueball) 7) < pip ~ nn \ n<f> ^ fo (quarkonium) 7) (2.3) 

< ^ r (^ fo (quarkonium) 7) (2.4) 

Thus if f (975) is a glueball, this branching ratio should be more than an order 
of magnitude smaller than it would be to a 0-like quarkonium. 
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If the /o is a quarkonium consisting only of nonstrange flavours, with a® its 
isovector quarkonium partner, these states will be OZI decoupled in the radiative 
decay. The OZI-violating production rate via a KK loop, viz. <fi — > ^KK — > 700, 
may be calculated. This calculation reveals some interesting points of principle 
which shed light on the role of finite hadron size in such loop calculations; this 
calculation will be discussed in the next section. 

Interesting questions arise in the case of qqqq or KK bound states ( "molecules" ) . 
The quark contents of these two systems are identical but their dynamical struc- 
tures differ radically. The situation here has its analog in the case of the deuteron 
which contains six quarks but is not a "true" six-quark bound state. The essen- 
tial feature is whether the multiquark system is confined within a hadronic system 
with a radius of order (Aqcd)^ 1 or is two identifiable colour singlets spread over a 
region significantly greater than this (with radius of order (fJ,E) 5 associated with 
the interhadron binding energy E for a system of reduced mass //). In the former 
case the branching ratio may be as large as 1CT 4 (see Ref. [5] and section 4); 
the branching ratio for a diffuse KK molecular system can be much smaller, as 
discussed below. 

The ratio of branching ratios is also interesting. The ratio of T(0 — > ^a ) /T(<fi — > 
7/0) is approximately zero if they are quarkonia (the fo being ss and the ao being 
OZI decoupled), it is approximately unity if they are KK systems, while for q 2 q 2 
the ratio is sensitively dependent on the internal structure of the states. This 
sensitivity in qqqq arises because — > Sj is an El transition whose matrix ele- 
ment, being proportional to Ee^, probes the electric charges of the constituents 
weighted by their vector distance from the overall centre of mass of the system. 
Thus, although the absolute transition rate for S = qqqq depends on unknown 
dynamics, the ratio of ao to fo production will be sensitive to the internal spa- 
tial structure of the scalar mesons through the relative phases in I — and 1 
wavefunctions and the relative spatial distributions of quarks and antiquarks. 

For example, suppose that the state's constituents are distributed about the 
centre of mass with the structure (qs)(qs), where q denotes u or d, and (ab) 
represents a spherically symmetric cluster. Then 




(2.5) 



and the El matrix element will be 



M ~ 



[( 



e u + ei) ± (e d + e 5 )] = e K + ± e K o 
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and hence the ratio F(<f> — > 7/o)/r(0 — ■> 700) will be unity. The quarks are 
distributed as in a fT-K^ molecular system (which is a specific example of this 
configuration) and only the absolute branching ratio will distinguish q 2 q 2 from 
KK. 

If the distribution is (qq){ss) then the matrix element 
M ~ [(e, + e ? -) - (e s + e g )] = . 

Here the quark distributions mimic 7T i] (in the do) or 7777 (in the /o). In this 
case the absolute branching ratios will be suppressed. Most interesting is the case 
where S = DD, where D denotes a diquark, i.e. where 

= -^=[(us)(us)±(ds)(ds)} (2.6) 




in which case 



so that 



M ~ [(e u + e s ) ±(e d + e 
r(0^ 7 a o ) 4 + 2,2 



= 9 



r(0^ 7 /o) v i-2' 
The absolute rate in this case depends on an unknown overlap between KK and 
the diquark structure; nonetheless the dominance of do over /o would be rather 
distinctive. For convenience these possibilities are summarised in Table 2. 



3 The KK Loop Contribution to — > S7 



The and the S (where S = a or / ) each couple strongly to KK, with the 
couplings g^, and g for <pK + K~ and SK + K~ being related to the widths by 

r(0 -> x+x-) 

and 

T{S ^ K + K-) 

for kinematical conditions where the decay is allowed. Hence, independent of 
the dynamical nature of the S, there is an amplitude M((f> — > S7) for the decay 
<fi — > S7 to proceed through the charged if loop (fig. 1), — > K + K~ — > 5(£) + 7 



K-4<+) 3/2 (3.1) 



" (3.2) 



167rm| 
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where the K ± are real or virtual and S is the scalar meson with four momentum 
£. The amplitude describing the decay can be written 

M(0(p, e0 )^5(£)+ 7 (g,e 7 )) = ^|-/(a,6)[(p-g)(6 7 - e0 )-(p- e7 )(g-6^] (3.3) 

where e 7 and (q and p) denote 7 and (f> polarisations (momenta). 

The quantities a,b are defined as a = —£-,b = so that a — b = ^2 j s 
proportional to the photon energy, and I(a, b) which arises from the loop integral 



is 



Ha, b) = - J— , - -f-^ifih ~ /(-)} + T^wMt) ~ 9(~)} (3-4) 
2{a — b) (a — b) z b a (a — b) A b a 



where 



m 

g(x) 



-(arcsin(2^)) 2 x > \ 
i[ln(^)-m] 2 x<\ 

(Ax - 1) 1/2 arcsin(^) x > \ 
I(l-4x)V2[l n (2±)_ X i7r ] x< \ 



V± = — (\±{l-Axfl 2 ) (3.5) 

Note that £ 2 may in general be virtual, though we shall here concentrate on the 
real resonance production where £ 2 = m 2 s with m s — 975 or 980 MeV. 

Even though Refs. [1-4] use essentially the same values for the couplings and 
other parameters, they obtain different results. Our results confirm those of Ref. 
[1] apart from a minor numerical error. Ref. [5] claims that the value of the loop 
calculation depends on the dynamical nature of the S. Since the coupling S — > KK 
is input from data it is somewhat surprising that the result can discriminate 
amongst models of the S. We confirm the numerical result of Ref. [5] and discuss 
its physical significance below. 

The resonant contributions to the <fi — > K°K 0/ -f branching fraction give a dif- 
ferential decay width 

dT |/(o,6)|W 



dk 2 " 4m A K 7i 4 



-X (3-6) 



where x is given by 



X \2M 2 m% {£ 2 - m 2 ) 2 + m 2 T 2 1 ' 
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Here £ 2 is the invariant mass squared of the final K°K° system, and hence the 
resonance. 

The limitations and problems in the existing literature concerning attempts 
to calculate the above are discussed in Ref. [11]. Here we shall briefly review 
the loop calculation in order to assess the existing literature and to highlight the 
novel features of the case where the S is a KK bound state with a finite size. 

Calculation of the integral I(a, b) 

Upon making the and K interactions gauge invariant, one finds for charged 
kaons 

H int = (e^ + g^f - 2eg^<j>^K (3.8) 

where A^, 0^ and K are the photon, phi and charged kaon fields, = iK\d^ — 
d )K. If the coupling of the kaons to the scalar meson is assumed to be simply 
the point-like one SK + K~ , then gauge invariance generates no extra diagram and 
the resulting diagrams are in figs. (1). Immediately one notes a problem: the 
contact diagram fig. (la) diverges. The trick has been to calculate the finite sum 
of figs, (lb) and (lc) and then, by appealing to gauge invariance, to extract the 
correct finite part of fig. (la). This is done either by 

a) (Refs. [1-3]) Fig. (la) contributes to A v 4>^g fJiU whereas figs, (lb) and (lc) 
contribute both to this and to p u q^A p '0^ . Therefore one need calculate only the 
latter diagrams, since the finite coefficient of the p v q^ term determines the result 
by gauge invariance. 

b) (Ref. [5]) These authors compute the imaginary part of the amplitude 
(which arises only from figs, (lb) and (lc)) and write a subtracted dispersion rela- 
tion, with the subtraction constrained by gauge invariance. This is also sufficient 
to determine the amplitude. 

In section 4 we shall consider the case where the scalar meson is an extended 
object, in particular a KK bound state. The SKK vertex in this case involves a 
momentum-dependent form factor f(k), where k is the kaon, or loop, momentum 
which will be scaled in f(k) by ko, the mean momentum in the bound state 
wavefunction or, in effect, the inverse size of the system. In the limit where 
R — > (or k — > oo) we recover the formal results of approaches (a, b) above, as 
we must, but our approach offers new insight into the physical processes at work. 
In particular, in this more physical case there is a further diagram (fig. (2d)) 
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proportional to f'(k) since minimal substitution yields 

f(\k-eA\)-f(\k\) = -eA.k^ (3.9) 

As we shall see, this exactly cancels the contribution from the seagull diagram fig. 
(2a) in the limit where q 1 — > 0, and gives an expression for the finite amplitude 
which is explicitly in the form of a difference M(q) — M(q = 0). This makes 
contact with the subtracted dispersion relation approach of Ref. [5]. 

First let us briefly summarise the calculation of the Feynman amplitude in the 
standard point-like field theory approach, as it has caused some problems in Refs. 
[2,3]. If we denote = [p u q^ - (p.qjg^H^^ms, q) (see eq. (3.3)), then the 
tensor for fig. (3) may be written (compare with eqs. 8 and 6 of Refs. [2] and [3], 
respectively) 

M f d A k (2fc-p)„(2fc-g)„ 

V em J {2nY{k 2 -m 2 K )[{k-q) 2 -m 2 K \[{k-p) 2 -m 2 K \ { ' 

We will read off the coefficient of p^q^ after combining the denominators by the 
standard Feynman trick so that 



(2tt)Wo Jo V Loo [(k-qy-pz) 2 -c + ie\ 3 1 ' 

where c = m 2 K — z(l — z)m 2 ^ — zy(m 2 s — m|). The p u q^ term appears when we 
make the shift k — > k + qy + pz to obtain 

h = n Jo dz si~ z dy yz[m2i< ■ z[i - z)m * ■ zy{m " ~ m * )] ~ i - (3 - i2) 

Note that m| < m 2 ^ and so one has to take care when performing the y integration. 
One obtains (recall a = m^/rri^, b = m 2 s /m 2 K ) 

H = l —{ t d Az(l -z)- (1 ~ Z{1 - Z)a) Inf 1 "** 1 "*^ ! 

AnHmlia-byio z [[ ] (a-b) [ l-z(l-z)a )l 

(a - b) Ji/ V+ z 

where rj± = |(1 ± p) with p = yl — A/a. (In performing the integrals, one must 
take care to note that a > 4 whereas b < 4 (which causes p 2 > 0,p% < 0)). Our 
calculation has so far only taken into account the diagram where the K + emits 
the 7; the contribution for the K~ is identical, so the total amplitude is double 
that of eq. (3.13) and therefore in quantitative agreement with eqs. (3) and (4) 
of Ref. [1]. Straightforward algebra confirms that this agrees with eqs. (9-11) of 
Ref. [5]. 



10 



Numerical evaluation, using m(fo) = 975 MeV and g 2 /Att = 0.6 GeV 2 leads to 

r(0 -> /o7) = 6 x KT 4 MeV (3.14) 

somewhat at variance with the value of 8.5 x 10" 4 MeV quoted in Ref. [1] fl. 
Ref. [5] does not directly quote a rate for (f> — > / 7- Instead, it quotes values for 
(f) — > 7/0 ~~ ^ 7 7r7r (for example) and claim that this depends upon the qq or q 2 q 2 
structure of the Jq. However, the differences in rate (which vary by an order of 
magnitude between qq and q 2 q 2 models) arise because different magnitudes for the 
fKK couplings have been used in the two cases. In the q 2 q 2 model a value for 
g 2 (fKK) was used identical to ours and, if one assumes a unit branching ratio 
for /o — > 7T7T, the rate is consistent with our eq. (3.14) (Ref. [5] has integrated 
over the resonance). In the case of the do, Ref. [5] notes that in the q 2 q 2 model 
the relation between g 2 (aoKK) and g 2 (ao7rr)) implies T(a — * ivrf) ~275 MeV. 
In the qq model, in contrast, Ref. [5] uses as input the experimental value of 
r(cio — > 7r?7) ~ 55 MeV which implies a reduced value for g 2 (ao7rrj) and, therefore, 
for g 2 (a Q KK): the predicted rate for (p — > 7<2 — * ■ynr] is correspondingly reduced. 

Thus we believe that the apparent structure-dependence of the reaction — > 
1S7 claimed in Ref. [5] is suspect. The calculation has assumed a point-like scalar 
field which couples to point-like kaons with a strength that can be extracted from 
experiment. The computation of a rate for <fi — > KK — > 76* will depend upon this 
strength and cannot of itself discriminate among models for the internal structure 
of the S. 

We shall now consider the production of an extended scalar meson [11] which 
is treated as a KK system (based on the picture developed in Refs. [15]). 

4 KK loop production of an extended scalar meson 

Suppose that K + and K~ with three momenta ±k produce an extended scalar 
meson in its rest frame. The interaction Hamiltonian H = gcf)(\k\)SK + K~ is in 
general a function of momentum. Now make the replacement k — > k — eA, and 

— * — # 

expand (f>(\k — eA\) to leading order in e; one then finds a new electromagnetic 
contribution 

H K +K- fol = -eg<p\k)k ■ A . (4.1) 
* However, J. Pestieau, private communication, confirms our value. 



11 



The finite range of the interaction, which is controlled by (f>(k), implies that the 
currents flow over a finite distance during the KK S transition: this current is 
the "interaction current" . The above current given by minimal substitution is not 
unique, in the sense that the transverse part e* 7 • j cannot be determined by the 
requirement of gauge invariance alone. However, it should describe the process 
under consideration accurately since the radiated photon is soft: the details of 
the interaction current are not important in the soft photon regime [25]. The 
effect of this form factor is readily seen in time ordered perturbation theory. (In 
this section we will work in the non-relativistic approximation. This suffices both 
to make our point of principle and to provide numerically accurate estimates for 
nonrelativistic KK bound states such as the f and a in the Ref. [15] picture. In 
general there are further time orderings whose sum gives the relativistic theory; 
see below.) 

There are four contributions: (#1,4 correspond to figs. (2a) and (2d), while 
if 2 ,3 correspond to figs. (2b) and (2c), where the 7 is emitted from the K + or K~ 
leg). We write these (for momentum routing see fig. (3)) 



#2,3 



#1 

#4 



2egg, Uk m ^ M ^f^ (4.2) 

yy *J D{E)D 1 D(q±) V ; 

2e m jA^ft (4.3) 

2em J d * k m±*M (4.4) 



where 



D 1 = m^-q-D{E) 

D(q±) ee m 4> -2E ± (4.5) 

D(0) ee m (j> -2E(k) (4.6) 

D(E) ee E + + E~ (4.7) 

and where E ± = E(k ± q/2) with E(P) the energy of a kaon with momentum 
P. Note that H 1 is the (form-factor-modified) contact diagram and H 4 is the new 
contribution arising from the extended SKK vertex. 

After some manipulations their sum can be written 

TT n ^ r rtl i(fc) f1 P-ik.qf, 1 1 <j)'(k)\k\ 1 , AnS 

H = 2e99 ^/ d 3 M ^ {1+ _^ ( _ + _ )} + ^JUi , (4 , 8) 
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If lim fc 2^ 00 (A; 2 (/)(/c)) — > | we may integrate the final term in eq. (4.8) by parts 
and obtain for it 

H 4 = 2egg&.$ f d 3 k^§-{-l - *°~ ( *£ } (4.9) 



D(0) E(k)D(0) 

This is identical to the q — > limit of H1+H2+H3, and hence we see explicitly that 
the g^ u term (i.e., the term proportional to e 7 .e^ as calculated above) is effectively 
subtracted at q = due to the partial integration of the (j)'(k) contribution, H4. 

If one has a model for 4>(k) one can perform the integrals numerically. For the 
KK molecule, the wavef unction 

M r ) = 1 u ^ (4.10) 

is a solution of the Schrodinger equation 

1 d 2 

{ — - + v(r)}u(r) = Eu(r). (4.11) 

rriK dr 2 

One may approximate (see Ref. [23]) 

v ( r ) = -440(MeV)exp(--(— ) 2 ) (4.12) 

2 r 

with ro = 0.57 fm. Equation (4.11) may be solved numerically, giving E = 
— 10 MeV and a ip(r) which for analytic purposes may, as we shall see, be well- 
approximated by 

^(r) = (^) 1/2 exp(-/ir); fj, = (4.13) 

where R kR ~ 1.2 fm (thus ^(0) = 3 x 10~ 2 GeV 3/2 ; see also Ref. [23]). The 
momentum space wave function that is used in our computation (see fig. (4)) is 
thus taken to have 

M = ^ . (4.14) 
0(0) (F + ^ 2 ) 2 1 ; 

The rate for T(0 — > S'y) is shown as a function of R^R m fig- (5)- The nonrela- 
tivistic approximation eqs. (4.2-4.9) is valid for Rrr ~ 0.3/m which is applicable 
to the KK molecule: for Rjck ~^ the fully relativistic formalism is required and 
has been included in the curve displayed in fig. 5. As Rkk ~^ and </>(k) -> 1 we 
recover the numerical result of the point-like field theory, whereas for the specific 

§ Actually, when k — > 00 the relativistic expressions of the next subsection are needed. These 
show that 4>(k) need only vanish logarithmically to obtain convergence. 
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KK molecule wavefunction above one predicts a branching ratio of some 4 x 10~ 5 
(width ~ 10~ 4 MeV). This is only | of the point-like field theory result but is 
larger than that expected for the production rate of an ss scalar meson (see Tables 
1 and 2). 



Connection with Relativistic Field Theory 



The nonrelativistic formalism is sufficient for describing the radiation from a 
KK molecule. However, it does not have the proper limit as Rkk — » 0; in this 
limit relativistic KK pairs are important in the loop integral. In this section we 
show how the relativistic formalism can be obtained from time-ordered perturba- 
tion theory and make contact with the relativistic field theory formalism of section 
3. The matrix elements for the time-orderings of fig. (6) are 

/d^k 1 
m K 1 ~ 2E + 2E_{E 3 + E + + E_) 

+ 2E + 2E„{m (f> -q-E+-EJ) ^ 

where the first (second) term corresponds to fig. (6a) (fig. (6b)) and E± is defined 
by E± = E(k ± q/2). Using E s = m 4> - q, 



m - q + E + + (m - q + E^) 2 - E% m <t> -q + E_ - E. 

and 

1 2E 1 

= +tzz : ~ v, F2 + ~ : F i F > ( 4 - 17 ) 



rritf, - q - E + - E_ (m^ - q - E + ) 2 - El - q - E + + E 

we obtain 

p r 



/d k -> 



+ ^t7tz^ — i 5 ~, m 1- ( 4 - 18 ) 



(2tt) 3 rvi '"" ~* L 2£_[(m - q + E_) 2 -E 2 

1 

2£+[(m^-(P£ + ) 2 -£ 

Analogously, M^ 4 , M^, and M4 are 

W) 

1 

2E + [( g + E + )2-^][(m^ + E + )2-E2] 



M£ = M» = +ieg g<p J ±^ <f>{\k\) 2e» [P - (k • qf] 
x [ 
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+ 2E.[(q - E„f - El][( m<t> - q + E„f - E\\ 

+ 2E + [( m<p - E + Y - El][(m^ - q - E + f - El] ] (419) 
d 3 k Ijfcl _ „ r 1 



+ 2B„[( m< - B )2 - £2] 1 (4 ' 20) 



where i?o is defined by £0 = E(k). 

In this way, we obtain "relativistic" expressions for the radiative meson 
decays. Matrix elements for the process a— d in fig. (2) may thus be written 

«f - c(t - 9/2 )S + , /2 - P ) < 4 - 21 > 

/■ d 4 /c e (j> -(2k + q-p) (2kY , inn . 

+ q/2)D{k - q/2)D{k + q/2 - p) v ' ' 

f dAk Ea,- (2k- q + p) (2kY . A nn . 



«r=W^^^ife^ (4 - 24) 

where D(k) is defined by 

D(k) = k 2 -m 2 K (4.25) 

and k = (0,k/\k\). In the particular case where 0(|fc|) = 1 and = 0, these 

reproduce the familiar field theory expressions of Refs. [1-5] and section 3. It is 
interesting to note the role that <fr'(\k\) plays in regularising the infinite integral. 

Define the matrix elements Mj (j = 1 — 4) by Mj = e 7 • Mj/ [iee^ ■ e^} and the 
decay width is then calculated by 

T(<f> -> S 7 ) = \M\ 2 , M = M 1 + M 2 + M z + M A (4.26) 

which reproduces the expressions in Refs. [1-5] and provides a check on our 
formalism. Eqs. (4.21-4.24), when evaluated numerically, give the decay widths 
shown in Fig. 5. In the limit R K R our numerical results agree with eq. (3.14) 
which was obtained by using the point-like field theory. 
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5 A Comment on the OZI Rule 



The calculations presented in this paper may have a bearing on one of the 
least understood characteristics of the low energy strong interactions: the Okubo- 
Zweig-Iizuka (OZI) rule [20]. If the a$ were a ^(mm — dd) state, its production 
in — > ao7 would vanish in "lowest order" in the quark model, with the KK 
loop contribution presumed to provide a small correction since such processes are 
OZI- violating {e.g., u — <fi mixing could also occur via such loops). We have seen 
that in the point-like approximation <fi — > ao7 would proceed with a branching 
ratio of order 1CT 4 via this loop process, as would /o = (uu + dd). If f = ss, 
a similar rate would be obtained from the KK loop, but now there would be a 
direct term which is supposed to be dominant. It is, however, easy to discover 
that this direct process would only produce a branching ratio of the order of 1CT 5 
(see Table 1). 

Our calculation provides some insight into this conundrum. If the KK system 
is diffuse, R K k ~ 2/m, then the loop calculation gives a branching ratio < 1CT 5 
(see fig. (5)) and the empirical OZI rule is good. Physically, the rate is suppressed 
due to the poor spatial overlap between the KK system and the 0. The point-like 
field theory does not allow for this: superficially the loops have a large magnitude. 
The essential observation is that the point-like calculation does not take into 
account the confinement scale, even though it is clear from our results that the 
dynamics can depend on it rather critically. 

Now consider a <fi and assume that S is an (ss) scalar meson, confined in 
Aqq D ~ 1 fm and connected by an intermediate state with quark composition 
qqss. If this multiquark system is confined in a length scale < ^qcd — l/ m 
(i.e., it is a "genuine" q 2 q 2 state and separate identifiable kaons are not present), 
then the point-like field theory calculations, which contain no intrinsic length, are 
superficially at least roughly applicable. The <fi branching ratio via the KK 

part of this compact system is then elevated above the 10~ 5 barrier. However, if a 
pure KK intermediate state forms, then it must occupy > 2Aqq D . The amplitude 
for the or a S(ss) to fluctuate to this scale of size would be small and it is this 
supression that is at the root of the OZI rule in this process. 

We see from this reasoning that the contribution of diagrams which correspond 
at the quark level to qqss loops really contain two distinct contributions at the 
hadronic level. These are first of all the diffuse contributions which can arise from 
hadronic loops corresponding to nearby thresholds, in this case from KK. Then 
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there are "short distance" contributions where approximating the qqss system 
as a KK system is potentially very misleading: a realistic calculation of such 
contributions would at least have to include a very large set of hadronic loops. 
A step in this direction has recently been taken in Refs. [26]. These authors 
have considered the loop contributions to, e.g., & — mixing in the 3 Po quark 
pair creation model, and found that there is a systematic tendency for the sum 
of all hadronic loops to cancel. In fact, they show that (in their model) the 
incompleteness of the cancellation of OZI- violating hadronic loops is precisely due 
to nearby thresholds. 

6 Conclusions 



There is still much thought needed on the correct modelling of the KK or q 2 q 2 
scalar meson and the resulting rate for <fi — > S7: the present paper merely makes 
a start by clarifying the present literature, making the first predictions for the 
production of a KK molecule, and pointing out the utility of the ratio of branching 
ratios as a filter. However, these results in turn raise questions that merit further 
study. For example, there are interesting interference effects possible between 
the ao(I = 1) and fo(I = 0) states which have not been considered. These two 
nearly degenerate states lie so near to the KK thresholds that the mass difference 
between neutral and charged kaons is not negligible: for example, their widths 
straddle the K + K~ threshold but only barely cross the K°K° threshold (at least 
in the case of the relatively narrow fo). 

Although there is clearly much to be done, it is already clear that there may 
be unique opportunities for probing dynamics in — > Sj and investigating the 
nature of the scalar mesons below 1 GeV. Moreover, we can already conclude that 
the branching ratio of <p Si will be between 10~ 4 and 10" 5 depending on the 
dynamical nature of these scalars and so will generate nugatory^] background to 
studies of CP-violation at DA$NE or other 0-factories. 
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Figure Captions 

Figure 1. The contact (a) and loop radiation (b,c) contributions. 

Figure 2. As fig. 1 but with an extended scalar meson. Note the new diagram 
(d). 

Figure 3. Momentum routing. 

Figure 4. Comparison between the exact momentum space wavefunction <f)(k) 
(solid) and the approximation of eq.(4.14); k is the relative momentum 
of the K and K. 

Figure 5. Y{4> Sj) in MeV versus Rkk i n f m - 
Figure 6. The two time orderings of fig. 2(a). 
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